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Motivation and Objectives 


¢Motivation: The probability of collision, P.,, 
between two Earth-orbiting satellites can 
often but not always be approximated 
adequately using the “2D P,” formulation 


«Objective: Implement an improved method 
to estimate collision probabilities 
—Use Coppola’s analytical “3D P,” formulation* 


—Validate using well-studied test cases and Monte 
Carlo methods 
—Compare 2D and 3D P, for archived conjunctions 


— *V. T. Coppola (2012) “Including Velocity Uncertainty in the 
Omicron Probability of Collision Between Space Objects”, AAS 12-247. Poe malls 
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Monte Carlo P, Estimation 


¢Collision probabilities can be estimated using 
Monte Carlo simulations 


—Computationally intensive, especially for low 
probability events 


‘Alfano analyzes twelve conjunctions in detail 
using Monte Carlo simulations 


—Benchmark test cases that can be used for 
validation of the 3D P, software 


—Includes cases where the 2D P, method both 
succeeds and fails 


a *S. Alfano (2009) “Satellite Conjunction 
Omicron Monte Carlo Analysis”, AAS 09-233. Doyle Hall | 5 


Coppola’s 3D Pc Formulation 


Coppola’ provides an analytically-derived formulation 
to calculate P, and its time derivative 


t+T( QAP 
Pa) 4 ; Ja 1D time integral 
‘0 dt 
dP : i 
7 -= f I (F, t ) d*t 2D unit - sphere integral 
t 
47 


¢ These integrals must be calculated numerically 


- *V. T. Coppola (201 2a) “Including Velocity Uncertainty in the 
Omicron Probability of Collision Between Space Objects”, AAS 12-247. POnean Le 


Coppola’s 3D Pc Formulation 


¢ Coppola provides an analytically-derived formulation 
to calculate P, and its time derivative 


t+T( QAP 
Pa) 4 ; Ja 1D time integral 
‘0 dt 


2D unit - sphere integral 


Analyzing the probability rate* provides new insight 
into the time dependence of conjunction risks 


*K.J.DeMars, Y.Cheng, and M.K.Jah (2014) “Collision 


Omitron Probability with Gaussian Mixture Orbit Uncertainty”, J. 


Guidance Control and Dynamics, 37(3) 979-985, 2014. POE 


Coppola’s Conjunction Time Bounds 


¢-Coppola also provides estimates for the 
bounding times of a conjunction 
—These often bracket the nominal time of closest 
approach (TCA), but not always 


«These bounds only provide a first-cut 
approximation for the limits of the numerical 
integration over time 


—These limits sometimes need to be expanded to 
bracket sufficiently the time(s) when dP,/dt peaks 


Omitron ‘V. T. Coppola (2012) “Evaluating the Short Encounter Assumption 


of the Probability of Collision Formula”, AAS 12-248. POye aie 


CARA’s Current 3D Pc Implementation 


¢ The 3D P, method is general enough to use 
—Gaussian Mixture Model state distributions 
—Complex dynamical motion models 
—Full 6x6 time-dependent state covariances 


¢ CARA’s current implementation uses 
—Single Gaussian ECI state distributions 
—The Keplerian two-body motion model 


—Full 6x6 ECl-state covariances, propagated using an 
analytically-derived state transition matrix* 


¢ Future plans include more advanced approaches 


"J. T. Horwood, et a/. (2011) “Gaussian Sum Filters for Space Surveillance: Theory and 
Simulations,” J. of Guidance, Control, and Dynamics, vol.34, p.1839—1851. 


= ‘e ‘ ‘ ans ioe A 
Omitron * S.W. Shepperd (1985) “Universal Keplerian State Tranisition Matrix,” Celestial Doyle Hall | 9 
Mechanics, vol.35, p.129-144. 


Illustration based on 
Alfano’s* test case #2 


NOTE: Actual 10 
surfaces are much 


larger and thinner 
y (m) “20 -20 x (m) 
2D P, assumptions: 3D P, assumptions: 
- Presumes straight trajectory . Trajectories are curvilinear 
(green) (black) 
- Presumes static covariances * Covariances vary throughout 
(blue) the encounter (pink, orange) 


= *S. Alfano (2009) “Satellite Conjunction 
Omicron Monte Carlo Analysis”, AAS 09-233. Doyle Hall | 10 


Approximations Used for 2D P, Estimation 


¢In terms of the relative position/velocity state vector and 
the associated 6x6 covariance matrix: 


ea _ +(t-t,,)V(t,, : 

x(t) = a 

V(t) V(t.) 

P= a . | A(t) Os, | 
Bit) C(t) 03, 03, 


¢Here t,, = TCA = the time of closest approach 
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Step-by-Step Relaxation of 
the 2D Pc Approximations 


Using CARA’s current 3D P, software, the three 
approximations used in the 2D P, method can be 
relaxed in a step-by-step manner: 


Coppola 1: Linear motion, A=A(TCA), B=C=0 
Coppola 2: Kep2Body, A=A(TCA), B=C=0 


= Coppola 3: Kep2Body, A=A(t), B=C=0 
Coppola 4: Kep2Body, P=P(t) 


Step “Coppola 1” employs all of the 2D Pc assumptions 

Step “Coppola 2” introduces Keplerian 2-body motion 

Step “Coppola 3” introduces time-varying position covariances 
Step “Coppola 4” introduces position+velocity covariances 
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Outline 


¢ Analysis of well-studied conjunctions 
Benchmark test cases analyzed by S. Alfano* 


= *S. Alfano (2009) “Satellite Conjunction 
Omicron Monte Carlo Analysis”, AAS 09-233. Doyle Hall | 13 


Validation using Alfano’s* Benchmarks 


¢ Benchmark case #3: 
—“Linear” case where 2D P, is known to be accurate 


—The 3D P, software correctly reproduces the 2D P, 
approximation, and Alfano’s benchmark P, value 


«Benchmark case #10: 
—“Nonlinear” case where 2D P, is known to be inaccurate 


—The 3D P, software correctly reproduces Alfano’s 
benchmark P, value 


«Other benchmark cases also analyzed (but not shown) 


= *S. Alfano (2009) “Satellite Conjunction 
Omicron Monte Carlo Analysis”, AAS 09-233. Doyle Halll| 14 


Alfano’s’ “Linear” Test Case #3 


=== Coppola bounds for y = 1e-16 m= Foster Pc = 0.10035 
Coppola (linear motion) Coppola Pc = 0.10034 
* MC (linear motion, N = 3e+07) * MC Pc = 0.100352 + 0.000058 
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These plots validate that the 3D P, software correctly reproduces both the 
Monte Carlo and 2D P, estimates, when using the 2D P, approximations 


= *S. Alfano (2009) “Satellite Conjunction 
Omicron Monte Carlo Analysis”, AAS 09-233. Doyle Hall | 15 


Alfano’s’ “Linear” Test Case #3 


=== Coppola bounds for y = 1e-16 === Foster 2D: Pc = 0.10035 
Coppola 1: Linear motion, A=A(TCA), B=C=0 Coppola 1: Pc = 0.10034 
lanl Coppola 2: Kep2Body, A=A(TCA), B=C=0 -===Coppola 2: Pc = 0.10034 
ae Coppola 3: Kep2Body, A=A(t), B=C=0 ='=="Coppola 3: Pc = 0.10034 
acon Coppola 4: Kep2Body, P=P(t) “enn Coppola 4: Pe = 0.10034 
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These plots validate that the 3D P, software correctly reproduces the 2D P, 
estimate, even when the 2D P, approximations are fully relaxed 


= *S. Alfano (2009) “Satellite Conjunction 
Omicron Monte Carlo Analysis”, AAS 09-233. Doyle Hall | 16 


Alfano’s “Nonlinear” Test Case #10 


mer= Coppola bounds for y = 1e-16 meee= Foster 2D: Pc = 0.29016 

Coppola 1: Linear motion, A=A(TCA), B=C=0 Coppola 1: Pc = 0.29016 
= Coppola 2: Kep2Body, A=A(TCA), B=C=0 -=-==Coppola 2: Pc = 0.28998 
lai Coppola 3: Kep2Body, A=A(t), B=C=0 ===" Coppola 3: Pc = 0.28108 
conan Coppola 4: Kep2Body, P=P(t) “ns Coppola 4: Pc = 0.36406 


dPc / dt (s”) 
= oi NO 
Cumulative Pc 


oS 
on 


These plots validate that the 3D P, software correctly yields different 


results as the 2D P, approximations are relaxed in a step-by-step fashion 


= *S. Alfano (2009) “Satellite Conjunction 
Omicron Monte Carlo Analysis”, AAS 09-233. Doyle Hall | 17 


Alfano’s “Nonlinear” Test Case #10 


=== Coppola bounds for y = 1e-16 === Foster Pc = 0.29016 
* Monte Carlo (N = 3000000) * MC Pc = 0.36380 + 0.00035 


Honan Coppola 4: Kep2Body, P=P(t) vv Coppola Pe = 0.36406 
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These plots validate that the 3D P, software correctly reproduces the 
Monte Carlo simulation, and that the dP,/adt profile has two blended peaks 


*S. Alfano (2009) “Satellite Conjunction 
Doyle Hall | 18 


= 
Omiczon Monte Carlo Analysis”, AAS 09-233. 


Alfano’s “Nonlinear” Test Case #10 


7a w= Foster Pc = 0.29016 
Alfano’s* original plot for 
g p ¢* MC Pc = 0.36380 + 0.00035 


cumulative probability Haven Coppola Pc = 0.36406 
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These plots validate that the 3D P, software correctly reproduces Alfano’s 
benchmark Monte Carlo results 


ial *S. Alfano (2009) “Satellite Conjunction 
Omicron Monte Carlo Analysis”, AAS 09-233. Doyle Hall | 19 


Outline 


«Motivation and objectives 

¢Overview of collision probability theory 

¢ Analysis of well-studied conjunctions 
=> ¢ Analysis of archived conjunctions 


2D vs. 3D results, repeating events, small-P, screening 
«Conclusions 
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Analysis of Archived Conjunctions 


¢-The 3D P, method has been applied to 80,453 
archived conjunctions 


—Actual events that occurred between 2016 April 1 
and 2016 June 1 


-Relatively few have appreciable 3D P, values 
—Only 11,211 (14%) have P, > 10°" 
—Only 5,761 (7.2%) have P,> 10° 
—Only 2,674 (3.3%) have P, > 10° 
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Analysis of Archived Conjunctions 


¢-The 3D P, method has been applied to 80,453 
archived conjunctions 


—Actual events that occurred between 2016 April 1 
and 2016 June 1 


-Relatively few have appreciable 3D P, values 
—Only 11,211 (14%) have P, > 10°" 
—Only 5,761 (7.2%) have P,> 10° 


—Only 2,674 (3.3%) have P, >= 10° 
This is the most important set for the CARA team 
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Archived Conjunctions with 3D P, > 10° 


¢For these, most 2D and 3D estimates were found to be 
relatively close to one another 
—71% have 2D P, and 3D P, within 10% of one another 
—85% have 2D P, and 3D P, within 30% of one another 


«But smaller subsets were found to differ significantly 
—5.6% have 2D P, and 3D P, separated by a factor of 3 or more 
—2.4% have 2D P, and 3D P, separated by a factor of 10 or more 


¢ The cases where 3D P, >> 2D P, are of significant 
concern to the CARA team 


—Threatening conjunctions could be overlooked when using the 
2D P, approximation 
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Repeating Events for Close Proximity Satellites 


¢ Objects persistently orbiting Archived conjunction involving two 
in close proximity can make satellites flying in close proximity 
repeated close approaches to means Concole boundetory= 1616 


one another Coppola 1: Linear motion, A=A(TCA), B=C=0 
—Satellites within formations or = -~"~" CORPO e Repeeeey enh eo 
clusters sienna Coppola 3: Kep2Body, A=A(t), B=C=0 
eee Coppola 4: Kep2Body, P=P(t) 


— These conjunctions can often 
be identified by their long 
durations 


— This can create multiple, 
blended peaks in dP,/dt 
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¢ These types of conjunctions 
explain some but not all of 
the archived cases that have 
3D P,>> 2D P, 


dPc / dt (s) x 10° 
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Isolated Conjunctions with 3D P, >> 2D P, 


mem Foster 2D: Pc = 1.32804e-05 


=== Coppola bounds for y = 1e-16 
Coppola 1: Pc = 1.32808e-05 


Coppola 1: Linear motion, A=A(TCA), B=C=0 


ae Coppola 2: Kep2Body, A=A(TCA), B=C=0 -===:Coppola 2: Pc = 1.32809e-05 
es Coppola 3: Kep2Body, A=A(t), B=C=0 ===" Coppola 3: Pc = 5.42875e-05 
a Coppola 4: Kep2Body, P=P(t) mn Coppola 4: Pc = 5.41930e-05 
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Archived conjunction where the 3D P, estimate exceeds 
the 2D P, estimate by a factor of about four. 
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Isolated Conjunctions with 3D P, >> 2D P, 


“=== Coppola bounds for y = 1e-16 w= Foster 2D: Pc = 3.64637e-15 
Coppola 1: Linear motion, A=A(TCA), B=C=0 Coppola 1: Pc = 3.64641e-15 
= Coppola 2: Kep2Body, A=A(TCA), B=C=0 -===Coppola 2: Pc = 3.60376e-15 
===" Coppola 3: Kep2Body, A=A(t), B=C=0 ===" Coppola 3: Pc = 1.17803e-05 
Sonia Coppola 4: Kep2Body, P=P(t) mn Coppola 4: Pe = 1.17956e-05 
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Archived conjunction where the 3D P, estimate exceeds 
the 2D P, estimate by several orders of magnitude. 
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A Screening Test for Small-P, Values 


° i j : About 80% of the archived conjunctions 
rsh anager age have (Mp)min > 10 and 3D P, < 3x10°1” 


Mahalanobis distances 
have small 3D Pc values 


¢ This correlation provides 
the basis for an efficient 
small-P, screening test 


¢ Applying this screening 
test eliminates the need 
to calculate 3D P, for 
~80% of all conjunctions 


Log, ,(3D P.) 


0 5 10 15 20 25 30 35 40 
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Outline 


¢ Motivation and objectives 
¢Overview of collision probability theory 
¢ Analysis of well-studied conjunctions 


¢ Analysis of archived conjunctions 


> ‘Conclusions 


= *S. Alfano (2009) “Satellite Conjunction 
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Conclusions 


«The CARA team has implemented Coppola’s 
3D P, formulation into software 
—Validated using Alfano’s benchmark test cases 
—Provides estimates for both P, and dP,/at 
—Provides insight into the time dependence of risk 


¢Archived conjunction analysis indicates that 


—Occasionally the 2D P, approximation can be very 
inaccurate 


—An efficient small-P, screening test can be used to 
speed processing for large numbers of conjunctions 
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Illustration based on Alfano’s test case #2 
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schematic Illustration of 2D P, Assumptions 


Actual nonlinear nominal Linear 
relative trajectory approximation 
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schematic Illustration of 2D P, Assumptions 


Relative position PDFs evolve in time 
NOTE: Actual 10 surfaces are much larger and thinner 
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Mahalanobis Distance 


The Mahalanobis Distance measures the difference 
between the positions of the primary and secondary 
objects, relative to the scale of their combined 
covariance: 


M(t) = (7 A r) 


where 
r=r(t)=r, -r, A=A(t)=A,+A, 


(relative position) (combined covariance) 
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Minimum Mahalanobis Distance 


¢ The Mahalanobis 
distance varies as a = Coppola bounds for y = 1e-16 
function of time during Coppola 1: Linear motion, A=A(TCA), B=C=0 
a conjunction amine Coppola 2: Kep2Body, A=A(TCA), B=C=0 
== == Coppola 3: Kep2Body, A=A(t), B=C=0 
i aa Coppola 4: Kep2Body, P=P(t) 


¢ The minimum value 
(Mp5) min Often occurs 
near the conjunction 
midpoint, but not 
always 


= 
on 


° (Mp) min Values vary 
significantly for 
different conjunction 
events 


on 


Mahalanobis Dist., Mp 
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0.6 0.8 1 
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Monte Carlo P, Estimation Procedure 


1. Sample the state PDFs for both the primary 
and secondary satellites 

2. Propagate the sampled states over the desired 
time span, checking if the separation becomes 
less than the combined hard-body radii 

3. If so, register a collision at the time the 
spheres defined by the hard-body radii make 
first first contact 

4. Repeat steps 1-3 to improve statistical 
estimation accuracy 


Omicron sonereiies 


